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ABSTRACT 


The development of a speed-of-light hard-kill weapon system for military 
applications represents a significant advancement in technology over present 
conventional kinetic weapon systems. Over the past two decades, the US Navy has 
successfully developed a megawatt-class chemical laser; however, under some maritime 
environments, the high power beam propagation was unable to delivery sufficient energy 
to kill a modern anti-ship missile (ASM) due to significant atmospheric absorption and 
the resulting thermal blooming process. A critical problem to resolve for the shipboard 
high-energy laser weapon systems is to develop a shipboard-compatible megawatt-class 
laser weapon at a wavelength where the atmospheric absorption is smallest. The 
megawatt-class Free Electron Laser (FEL) has significant advantages over conventional 
weapon systems and other chemical high-energy laser systems. Infinite magazine, rapid 
response, and wavelength tunability make the FEL a suitable and desirable shipboard 
weapon system. 

This thesis divides into four chapters. Chapters I and II introduce the FEL and 
background theory of the FEL. Chapter II explores the analysis of the LANL 
Regenerative MW FEL Amplifier design and optimizes its efficiency. Lastly, chapter IV 


summarizes the feasibility of achieving the desired efficiency. 





TABLE OF CONTENTS 


eNO ODAC FI COIS] Se nn oP ] 
ieee) Fe Eo eyes Da ebel ev) be tI) CO GAN ooo sso sackcsdaaucessesovtscovessosusbbeclsbupeyecevcteoeas: 1 

B. SHORT HisTOR OF PREESELEGTRON LASER .............008..4961...0000:...08%. 2 

Ce BASIC TREE BIBBOTRON WAGER THEORY sccsccistiscsseseesssosesssorssersteesssmcdvueeree S 
TERN GCIPL ES OG@ RE Res dilO NO Te gE cose cen seerayeertoeseteer se nn nov essenzvevs 3 
A BEECTRON TRAJECTORY IN THE WNDULATOR © ...c.c0 nee 5) 

B. FEL RESONANCE —- ELECTRON AND PHOTON RACE Lu... eee 6 

C.FPE ND UT UNA OUAMIMION: ....<2i3 ee eee cetera acetate 8 

DP O PAG le Woe te UW AION: eater eetereerts arto tss dents ee eee 10 

E. GAIN DEGRADATION DUE TO ELECTRON BEAM QUALITY 0.000. 15 

Bs PABPERE DP WINIDUEAT OR occas tei ea anton, eee Ws 

lil, SIMULATIONS OF THE REGENERATIVE MW FEL AMPLIFIER ....0 0. 21 
A». MW-CLASS FEL CALCOLA TIONS Wicccteccetieecsesteeisse sce enn | 

B. LANL. REGENERATIVE MW GEL AMPEIBIRR 2 ne eee 21 
.-1D Efficiency Optimization 2c ceesccescics ee 24 

2. 3D Diifmaction Optimization: .... rojas 52 

IW CONCLUSIONS o.coiicececccsecsetetecntsccce-oseueh oer 45 
1eTST, OF REPERENCES ow csssss., cccceiesscasescige: sotto. seek 47 
INITIAL: DISTRIBUTION LIST 3... tee 49 


Vil 





ACKNOWLEDGEMENT 


The author gratefully acknowledges the guidance, advice, and attention to detail 
of Professors W. B. Colson and R. L. Armstead. The author would also like to thank to 
Dr. R. K. Wong for his efforts in assisting computer simulations. Their patience, insight, 
and technical assistance were invaluable. 

Finally, this thesis would not have been possible without the love and support of 


my wonderful wife KimChi, and my beautiful daughter Alisa. Thank you! 





I. INTRODUCTION 


A. DEFENSE AT THE SPEED OF LIGHT 


In today's technology, weapon systems rely predominantly on kinetic kill, striking 
a missile with another missile or with several bullets. These weapons are relatively slow 
and require a longer detect-to-engage sequence when factoring in tracking time, launch 
time and time of flight to intercept. In searching for faster and faster weapons and ever- 
decreasing response times, Directed Energy Weapons (DEW) have been on the U.S. 
military's "love to have" list for several decades. DEWs are attractive for several reasons, 
and the most predominant reason is that they offer speed-of-light defense. 

In 1971, the U.S. Navy established the Navy High Energy Laser (HEL) Program 
with the intentions of developing speed-of-light hard-kill weapon system. Initial 
development focused on building a CO, gas dynamic laser technology, but a few years 
later its emphasis shifted to continuous-wave deuterium fluoride (DF) chemical lasers 
because they propagate far better in the maritime environments. In the early 1980s, the 
megawatt class Mid-Infrared Advanced Chemical Laser (MIRACL) was built by TRW 
for the Navy as the highest average power continuous wave (CW) laser in the west. The 
SeaLite Beam Director (SLBD) was built by Hughes Aircraft to accommodate and direct 
the MIRACL beam. The MIRACL and SLBD have successfully tested at White Sands 
Missile Range, New Mexico, in engaging and destroying several types of missiles in 
flight [1]. 

Although the MIRACLE/SLBD system was successfully tested against crossing 
targets in a fleet defense scenario, its lethality was not sufficient in some ship self-defense 
engagement scenarios due to insufficient crosswind [1]. Insufficient crosswind will cause 
the laser beam heating the air in its own path resulting in a defocusing optical lens known 
as thermal blooming [2]. The thermal blooming effect can severely reduce the High 
Energy Laser Weapon Systems' (HELWS) lethal range. As a result, it's necessary to 


replace the MIRACL with a megawatt-class laser device at a wavelength where 


atmospheric absorption is smallest. An ideal solution is the megawatt-class Free Electron 


Laser (FEL) that is tunable to the optimum propagation wavelength. 


B. SHORT HISTORY OF FREE ELECTRON LASER 


The first generation of coherent radiation from free electron beams was 
microwave tubes that relied on slow-wave structures, and their development received a 
significant impact from the radar development during World War II. In 1951, Motz 
showed that an electron beam propagating through an undulator magnet can be used to 
amplify radiation. After 1960, research on short wavelength lasers has been dominant, but 
they require an atomic or molecular medium to operate, resulting in a limitation on 
wavelength tunability. In 1971, J. M. J. Madey of Stanford University brought attention to 
what he proposed as the 'free electron laser’. Relying on research of synchrotron radiation 
sources, Madey conceived a device that would lase in the visible using a beam of 
relativistic electrons as the source. In 1976, the first successful free electron laser (FEL) 
experiment was demonstrated by Madey and his co-workers at Stanford University, in 


that they measured 7% gain from an FEL configured as an amplifier at 104m 


wavelength. One year later, another successful FEL operation of the same FEL configured 
as an oscillator. These successful experiments created a large interest in FEL research. 
Since then, scientist and experimental groups around the world built FELs that operate at 
frequencies ranging from microwave to the UV [11]. Three other FELs were built soon 
after. The first was built at Laboratoire pour !'Utilisation du Rayonment 
Electromagnetique in Orsay, France, which operated in the visible range. The second 
involved a team from TRW at Stanford using the Superconducting Accelerator (SCA), 
lasing in the near infrared (IR). The third was built at Los Alamos that lased in the mid- 


infrared [3]. Now, there are many FELs used for research around the world. 


c. BASIC FREE ELECTRON LASER THEORY 


The FEL is a device that transforms the kinetic energy of a relativistic electron 
beam into electromagnetic (EM) radiation [4]. This is accomplished by sending a 
relativistic electron beam from a linear accelerator, microtron or a storage ring etc., 
through an alternating magnetic field produced by a device known as an undulator as 
shown in Figure 1.1. The interaction of the negatively charged electrons with the spatially 
periodic magnetic field induces a transverse acceleration on the electrons that causes 
them to exchange energy with a copropagating radiation field. The kinetic energy 
extracted from the electron beam is converted to electromagnetic radiation. The energy of 
the electron beam and the period of the undulator magnet determine the wavelength of 
radiation. This system provides a powerful source of tunable, coherent, electromagnetic 
radiation that can operate in a wavelength range from about 1 cm to the X-ray region. 
This system can also provide continuous wave (CW) operation, or very short pulses in the 


nanosecond, picosecond, or even subpicosecond region. 


Electron Beam 






Undulator 


Resonator Mirror 


Figure 1.1. Schematic diagram of a free electron laser. 





I]. PRINCIPLES OF OPERATION OF FEL 


A. ELECTRON TRAJECTORY IN THE UNDULATOR 


For better understanding of the physical mechanism underlying the operation of an 
FEL, we first consider the interaction of a free electron in the undulator and an 
electromagnetic wave. The Lorentz force equation and the equation for the energy 


exchange between the electron and the electromagnetic wave are 


poe {es ud 8) 2.1) 
Se eae. < . 
dt C : 
d(ymc? 7 
Eo. B. OF) 


where e, m and ¥ are the charge, mass and velocity of the electron, respectively; y is the 


relativistic factor of the electron, yv=1/ v(1 3. F B=v/c; p is the relativistic 

momentum of the electron, p=ymv ; E is the electric field of the EM wave; B is the 

magnetic field ; ync* is the electron energy; and c is the speed of light in vacuum [5]. 
Considering an electron with velocity }= fc and a transverse electromagnetic 


wave propagating colinearly, in the z-direction in a vacuum, then equations (2.1) and 


(2.2) become 





dy B)_ es 
zs, =-— (E+ x8), (2.3) 
On iri 
5 alae (2.4) 


Assume the undulator has a helical magnetic field of the form B= B[cos(k,z),sin(k,z),0], 
where k,=27/A, is the undulator wavenumber, and 4, is the undulator wavelength. 
The corresponding circular-polarized, plane-wave optical field is then 

E,= E(cosy,—siny,0), and B. =E(siny,cosy,0), 
where y=kz-—at+@ . E is the electric and magnetic field amplitude in cgs units, ¢ 1s 
time, @ is the optical phase, k =@/c is the optical wavenumber, and @ 1s the optical 


angular frequency. Substituting the fields into Lorentz force equations yields 


5 


AED r eB cosy. siny.0) + B_ B(-sin(k,2),cos(k,2).0)}, (2.5) 


d 
eee (GO, cosy— BP, sul)" BPD sinGaa)— Prcos(ezyr (2.6) 
dy e 
and ——=--——E(£,cosy-£,smy), (ay) 
dt mc ‘ 


where 6 =6,+f., B, =f,i+B,j,and PB, = B.k. For relativistic electrons, B, ~1, 
hence E(1—- £8.) « f, B, so in the right-hand side of (2.5) the first term is ignored. The 
second term 1s 


d(yp , ye eBp. 


(—sin(k,z),cos(k,z),0). (2.8) 
dt mc 
Integrating the transverse equation (2.8) and assuming a perfect injection of electron 


beam into the undulator cavity, which makes the constant of integration zero, yields 
= eB 
yp == (=) (cos(k,z), sin(k,z),0). (2.9) 
mc 
Divide both side of (2.9) by y to get 
5 K 
Bp (2=- py cose 2) suit, 2):0)) (2.10) 


where the undulator parameter is K=eB/k,mc* =eBA,/2nmc’ . Equation (2.10) 


describes the transverse motion of electrons in the helical undulator [6]. 


B. FEL RESONANCE -- ELECTRON AND PHOTON RACE 


Accelerated electrons traveling at a relativistic speed along the longitudinal axis of 
the undulator of the FEL radiate optical energy [6]. In order for the FEL to have gain (1.e. 
the fractional change in optical power in a single pass through the undulator), a net 
transfer of energy from the electron beam to the optical beam must occur. In order for an 
optimum exchange of energy from an electron to the optical field to occur, the electron 
must oscillate in phase with the optical field. This optimum energy exchange occurs when 
the electron falls behind the optical fields by a distance of one optical wavelength within 


the span of one undulator period [6]. 


Speed of electron £,c 
-— 


Speed of light c 


h-- +4 
Light wavelength 


t=0 Undulator wavelength /, t=4,/ B,c 


Figure 2.1. Electron-Photon race. 


Consider an optical wave that travels at speed c and passes over an electron traveling at 


speed £.c as shown in Figure 2.1. At resonance condition, one wavelength of the optical 


wave passes over an electron while it travels through one undulator wavelength 4, , so 


ee ne ; 
that light wavelength A=(c— £,c)t. Insert t=2,/B,c yielding inl A, . Substitute 


Z 


2 


Bp} =— into y? =1- B* =1-B; - BY yields 
Y 


(l+ K*)y? =1-£?. (2.11) 





lake 


Sineess- =|, |p = pues 


a sail : (2212) 

oy 
The equation (2.12) is the resonance condition that demonstrates the tunability of the 
FEL. By adjusting the undulator parameter K, or the energy of the electron beam 
(y-1)mc* or the undulator wavelength /,, the wavelength of the optical field can be 


tuned to a desired wavelength within a prescribed bandwidth [7]. This is a significant 


advantage over conventional lasers that are typically constrained to a fixed wavelength. 


C&: PENDULUM EQUATION 


This section derives the equation of motion of the electrons resulting from the 


interaction with the optical field. Substitute (2.10) into (2.7) to get 
y SEC ds (2.13) 
ymc 


where 

C+p=pwthz=(kt+k,)z-at+¢. (2.14) 
For -z/2<€+6<a/2, the cosine function in (2.13) is positive and the electrons will 
gain energy from the optical field, whereas for 7/2<¢€+4+@<3z/2, the cosine is 
negative and the electrons will lose energy to the optical field. Differentiating (2.11) with 
respect to time ¢ gives 27” A +K° )= 26. B,. Relating y to B, gives 


2 , 
Yells (2.15) 
YAN hee) 
Differentiating (2.14) once yields ¢ =(k+k,)2-—w=(k+k,)cBh,-@, and 
differentiating twice yields ¢ =(k+k, cB... Insert the resulting expression for B. into 
(2.15) toreet 
ein. a (2.16) 
y (1+K’*)(k+k,)e 


Since 2 « 4, then k, « k, so (kK +k, )c = kc = ow. Using this approximation and the 


resonance condition (2.12), equation (2.16) becomes 


Y , 





= 2a 
y (1+K’)o 2a, G17) 
Ce ot eke | —_ 
So, > — =— => cos(¢ + @), then the electron equation of motion is 
20, YY yme 
elem 
¢ =— 7 cos(¢ + @). (2.18) 
vy “Mec 


Equation (2.18) describes the electron phase dynamics in the form of the pendulum 
equation [6]. 
Consider time evolution through the undulator L/ £,c for an average electron, 


where L is the length of the undulator. Introducing relevant dimensionless time 


t= f,ct/L =ct/L, we see that t will vary between 0 and | along the undulator length. 


Substituting dt = (L/c)dr into the pendulum equation (2.18) gives 
¢ =v=|alcos(f +¢), (2.19) 


where the notation (..) =d(..)/dt, 





a| = 4nNeKLE/y’ mc’ is the dimensionless optical 
field strength, and N = L/2, is the total number of undulator periods. Equation (2.19) is 


the electron pendulum equation that governs the phase-space motion of the electrons 


within the undulator and under the influence of the optical wave [6]. 


D. OPTICAL WAVE EQUATION 


An electromagnetic wave will propagate in the absence of a current density 


according to 


[v" Se a = lac. a (2.20) 
Com 


where the Laplacian is V* = oe 4 om + ole , F is the vector position, ¢ is time, and A is the 


optical vector potential [5]. The optical electric and magnetic fields can be expressed in 


terms of the vector potential by the following relationships 


E = aa (2.21) 
GG 
and B=VxA. (222) 


Assuming that the optical beam is a circularly polarized plane wave, then the vector 


potential has the form 


a t) 


4 





(sin y, cosy, 0)=Re {i Birnee«| ; (2.23) 


where @ (7,1) =E(7, tHe?’ , w =kz—at+¢ is the plane-wave phase, é =(—i,1,0) is the 
polarization vector, k = 2zc/q@ is the optical wave number, @ is the optical frequency, 


and @(7,f) is the optical phase [7]. For convenience, we write 
A= 8G tee" (2.24) 
where @ = kz — wt. Differentiating (2.24) with respect to ¢ gives 


a 2(¥ ios}, 025 


where (.)=0/..) /t . Differentiating (2.25) with respect to t gives 


oA 


a =<(8- 2iw€ -a sia (2.26) 


where es )=0° (..)/ Ot’ . Differentiate (2.24) with respect to z gives 
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a a a i (2.27) 
Oz ok 
where (..)’=0(..) / 0z . Differentiate (2.27) with respect to z gives 


6° A _ wt 1@ 

~~ <( +2ik B'-P Fh", (2.28) 
where (..)” =0°(..)/@z’ . Since the laser beam consists of a narrow range of frequencies, 
the variation in the electric field amplitude and phase with respect to space over an 


optical wavelength 4 is small, so that &’« k@ and ¢'«k@. The variation in the electric 
field amplitude and phase with respect to time over an optical period is also small, so that 


€«oE& and ¢«od [7]. Therefore, &” is negligible compare to k@’, and & is 


negligible compare to & , then (2.26) can be simplified to 











a _ | 210 & - 08), (2.29) 
and (2.28) can be simplified as 
2 
a a (2ik &- kB). (2.30) 
Substitute (2.29) and (2.30) into the left hand side of (2.20) yields 
l eos ee'* ale... 
[v" 45 lar. : vis sn(s tal (2.31) 
O° 2 
where ve" ee . Hence, 
vi sain S41 5) @ =0. (2.32) 
Oz C Of 


The first partial second derivatives with respect to x and y describe how the beam changes 
due to diffraction as it propagates. The last partial derivatives with respect to z and / 
describe the propagation of a plane wave along the z axis. 

Introduce dimensionless time t=ct/L so that rt varies from 0 to 1] along the 
undulator length. It's convenient to follow the wavefront with a new coordinate 


Z =z-—ct. With this coordinate change, 


ars( gy, (Be) 2s 2 233) 


a \ale \ajar - OZ 
since 02 /0z =1 and Or/oz =(. Also, 
|i goes 
Bn (F)2+(2)2 i ae tS (2.34) 
Cromer c| or) Oz \oryror G7 anor 
since 0z/0t=-c and Or/ot=c/L. Therefore, the propagation operator in (2.32) 
becomes 
eo | as can 
Oz semel) OZ G27 1NOr } Pree 
Hence, the full wave operator is 
viva S415) =| vi 42k S| (2.36) 
CZ CU INOT 


Introducing the dimensionless transverse variables ¥ =xVk/2L and ~=yVk/2L, so 








6° 6° a o° ) me ek || 
that —-=(k/2L)—~ and —~=(k/2L)—. Hence V7 =(k/2L) —-+—, |= 
ax ox? dy oy x Oy 


2 Z 
(k / 2L)V*. , where V+ <= - — . Substituting into (2.36) and (2.32) yields 
Pas +2 & =0. (2.37) 


Equation (2.36) is known as the parabolic wave equation [7]. In dimensionless 


parameters, the optical field is represented by a= lale’* where 
la| = 42 NeKLE (y*mc*) « E ,so & and a differ by a multiplicative constant. Therefore, 


exchange @ for a, the dimensionless parabolic wave equation (2.37) becomes 


= O 
= Vee ain Ot 2.38 
E i 2 ( 


In the presence of the electron beam current density, the optical wave equation 


(2.20) is 


ech ee An - 
iv -45 ae = =e (39) 


where J, is the transverse electron beam current density [8]. Substituting (2.30) into 


(2.39) yields 





= vison S412 pee ES (2.40) 
k Oz 


A 


Performing a dot product with é =(i,1,0) on both sides of (2.40) and using 


é* -6 =(i1,0)-(-i,1,0) = 2, yields 


Ip + if 442 | @ = -="),. Ca (2.41) 
2 Of Ok 


The total transverse electron beam current is the sum of all single particle currents. The 


single particle current for one electron in the beam is J,, =—ecf,5°(¥—7) [8] where 


th 


F is the position of the i” electron and 6“) (¥,—7) is the three-dimensional Dirac delta- 


function. The total beam current is 
Ju Jj See, ale (X 7, )- (2.42) 
From (2.10), the transverse motion is £, =(—K/y\Xcosk,z,sink,z,0)= 


way ° ° A —_ K oA —1k 2 
Re (- K/y)iée“* } For convenience, we write £, =——iée™“’ , so that 
y 


J, Be =iecK(e’ - ao ee 68-7). (2.43) 
We can simplify the sum » by assuming a local electron density p, so that (2.41) 


becomes 


gate 
ie Vit i SAL) & = ~2nzieKp( — (2.44) 
2 oz cot y 
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where ¢ =k,z+q@_ is the electron phase in the combined undulator and optical fields, 
and fees) is a Weighted average of sampled electrons at a specific site which replaces the 


summation over all electrons within one optical wavelength. Define the dimensionless 


beam current 7 =8N(ez KL)’ p/(y’ mc’) , the wave equation becomes 


= O _ “pe 
2 +S Jars, 2.2 )=— Jf aa a) 


In some cases where the optical wave diffraction is not significant, we can ignore V‘ in 


(2.45), then 


gee jhe | (2.46) 
The optical wave equation (2.46) couples the electrons to the copropagating optical wave 


by relating the electron phase to the change in the optical field. The equation (2.46) can 


be expressed in terms of amplitude and phase, 


lal = —j(cos(¢ +¢)), (2.47) 
° _ ji ; 
and g¢= ia sins + @)). (2.48) 
a 
Equation (2.47) reveals that when the electrons become bunched near (¢ + ¢) = 7, the 
optical field amplitude will grow, and when bunched near (¢ + ¢) =~ 0, the optical field 


amplitude will decrease. Equation (2.48) describes the optical phase evolution. 
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E. GAIN DEGRADATION DUE TO ELECTRON BEAM QUALITY 


In the design of an FEL, there is often a trade-off between high beam current and 


good beam quality as determined by the accelerator. The dimensionless current 
j«IN°A'* where J is the actual beam current , N is the number of undulator periods, 


and / is the optical wavelength. Usually the optical wavelength 4 is fixed for a given 
application so that / is maximized by adjusting current / and the undulator periods N. 
Increasing / tends to degrade the beam quality from the accelerator, whereas increasing N 
tends to increase the FEL sensitivity to the beam quality because of narrowing the gain 
spectrum bandwidth [9]. 

Emittance € describes one aspect of electron beam quality. A real electron beam 
has an emittance ¢ describing how electrons enter the undulator with a range of angles 


@. with respect to the undulator axis, and at a range of transverse positions x, slightly 
off axis. The electron beam may also have a small energy spread Ay/y. The energy 
spread and emittance of the electron beam translate into a spread in phase velocities in 
phase-space. The phase velocity spread due to energy spread is given by [10] 

Av =4zNAy/y, (2.49) 
and the phase velocity spread of an electron beam due to emittance is given by 


272 N 


Av =- 5 
I+K 





(K?k2x? +7767). (2.50) 


Electrons enter the undulator with a spread in phase velocities due to a spread in 
energies, and are randomly distributed in both the longitudinal and transverse directions. 
The typical distribution of electron phase velocities is assumed to have a Gaussian form 


about initial phase velocity v,, 


—q’? 122 
I6(9) fee ) 
° V2 O- 


where 0, =4aNAy/y_ is the standard deviation and q is the phase velocity away from 





(2.51) 


v, distributed [9]. The characteristic function F(t) of the distribution f,(q) is defined 


dS 
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Fo(t) = |fo(qve dg. 
Evaluating the above integral equation for a Gaussian yields 

F(t) mentor? 
If the beam quality is perfect so that o, = 0, then F(t) =1. However, when the beam is 
not perfect, the characteristic function Fe (z) decays over the length of the undulator in 
a characteristic time on the order of 1/a, [9]. 


Poor beam quality reduces the single-pass FEL gain by degrading electron 
bunching. As an example, consider an FEL design with N = 25 periods having a 


dimensionless current density of 7 =10 and initial phase velocity of v, =3 operating with 
an initial dimensionless optical field a, = 1. In Figure 2.2, electron beam quality is 
assumed to be nearly perfect, i.e. o, =1. Some electrons gain energy and move ahead of 


the average flow, while other electrons lose energy to the radiation field and move back 


behind the average flow causing visible spatial bunching. 





Figure 2.2. Phase space evolution for a nearly perfect electron beam quality. 
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In Figure 2.3, the longitudinal energy spread increases to o, = 6, and this 


degradation in beam quality causes much less bunching when comparing to Figure 2.2. 
The poor beam quality is responsible for a decrease in final gain to Gy = 0.10534 and a 


smaller optical phase shift @. 
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Figure 2.3. Phase space evolution for low gain with electron beam degradation. 


F. TAPERED UNDULATOR 


The characteristics of the FEL interaction can be altered by changing the 
undulator properties, i.e. undulator polarization, wavelength, or field strength, along its 
length [5]. As an FEL reaches saturation, the electrons have given up enough energy so 
that they are no longer in resonance with the optical wave, which decreases the gain. 
Tapering changes the undulator's resonance condition along its length so that the 
electrons remain in resonance longer, hence increasing the overall gain in strong fields 


[7]. From the resonance condition (2.12), in order to extract energy oul of the electron 


Ia 





beam, i.e. reducing y, and convert into radiation field and still maintain resonance 
condition, we can either taper the undulator wavelength /,, i.e. reducing 4,, or taper the 
undulator parameter K, 1.e. reducing the undulator field strength B. However, tapering K 
is eaSler since the undulator gap needs to be increased to decrease B. Tapering essentially 
has the same mathematical result as accelerating the electrons. This artificial acceleration 
is represented by a dimensionless parameter 6. A tapered undulator represents an 
advantage over the untapered undulator design when the artificial acceleration 6 exceeds 
the deceleration that can be obtained without taper in strong optical fields. When the 
undulator wavelength is decreased, the artificial acceleration 6 is given by 


0 =—20NAA,/A,. When the undulator field strength is decreased, 


5 ~-4nNK’AB/ B(1+ K*). With tapering is included in the undulator, the pendulum 


and wave equations become 
6 =v=6+l|a\cos(¢+¢), (2.52) 
and pans = 


Where ¢ = [ k(2)dt' + ke - at [10]. The wave equation remains the same, but the 


simple pendulum equation includes an additional constant torque due to the artificial 
acceleration 6. 


The necessary and desirable criteria for tapering 1s ld 0 > lal” = 27 ine 


condition on the left is for trapping electrons in the tapered phase space; the middle 
condition 1s that the tapered acceleration exceeds the natural untapered deceleration; and 
the condition on the right is that the tapered FEL works in the strong field regime [10]. 


Electrons initially near the phase ¢, ~0O are accelerated by both the strong optical 
field |a| and the artificial "torque". The electrons initially near ¢,~ undergo a 


cancellation between the optical field and the torque, leaving them trapped in closed 


orbits. The taper is effective because electrons near the phase for gain ¢ ~z are trapped, 


and continue to contribute to the interaction. Electrons near the phase for absorption 
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¢ x0 are taken away from resonance and eventually stop interacting. The imbalance 
leads to a net gain [10]. 

Compared to the untapered undulator, the tapered undulator is more efficient in 
strong optical fields, but has smaller gain in weak optical fields. The efficiency of a 
tapered undulator is estimated at 7, =d/(8aN). For example, if 6 ~ 2807 for an 
undulator with 23% taper in wavelength over N = 300 periods, then the efficiency is 


estimated at 7, ~70/(2.N) which increases by a factor of 70 over the natural efficiency 


n =~1/(2.N) [10]. Tapering is a design to obtain higher overall gain and efficiency in 
strong optical fields. 

An example of a tapered undulator in strong optical fields is shown in Figure 2.4 
with an untapered undulator shown in Figure 2.5. Both designs are based on an undulator 


length of L = 6 m corresponding to N = 300 periods. The dimensionless current density is 


J = 12000 » z in the high gain regime, and initial field a, =150 >»z in the strong optical 


field regime. In Figure 2.4, the tapering rate 6 = 2607 turns on at 7, = 0.19 along the 


undulator resulting in a final gain of Gre = 556 and FEL extraction efficiency of 


n #12.5%. Initially the gain increases exponentially until 7, =0.19 where the taper 


turns on and saturation is reached. With tapering, the gain continues to increase slowly 
instead of decreasing after saturation, as shown in the upper right-hand portion of the 
Figure 2.5. The phase velocity of the untrapped electrons near v ~740 1s shifted by the 
tapering acceleration 6 whereas the trapped electrons remain near resonance due to 
strong optical fields. The lower right-hand portion of the Figure 2.4 is the optical phase 
shift, and this phase shift Ag ~ 1.47 over the length of the undulator is a characteristic of 


high gain. 
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Figure 2.4. Phase space evolution for high gain in strong fields with a tapered undulator. 





Figure 2.5. Phase space evolution for an untapered undulator. 
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Il. SIMULATIONS OF THE REGENERATIVE 
MW FEL AMPLIFIER 


A. MW-CLASS FEL CALCULATIONS 


In order to destroy an incoming missile, the laser power required on the missile 
surface is about 10 kW/cm” over a w,* 5 cm radius spot for a two second duration. The 
extinction coefficient due to aerosols at sea level is @ ~ 0.03 km’, and through e“ 
describes the absorption loss at Alum wavelength over a distance z. Let P be the 
Transmitted Power that must leave the ship to destroy a missile at a range R, then the 


power at the target is P,=Pe™. If the range of the missile is 5 km, the laser power 


which must leave the ship is then P, = {10 kW/cm’] [x(5 cm)*] ce" *™ = 910 kW. 
As an example, consider an FEL with N = 300 periods using a beam of 100 MeV 
electrons with a peak current of 400 A, then the peak electron beam power P = 400 


Ax 100 MV ~ 40 GW. The micropulses generated by the FEL have a length of 0.001 m 
and are separated by 6 m, then the fraction D = 0.001/6 denotes the duty cycle of the FEL. 


The average power generated within the FEL 1s given by P=P xDx6.7 MW. The 
efficiency necessary to destroy an ASM is n=P,/P ~13%. Therefore, the single-pass 
FEL efficiency yn needed to supply the Laser Power to destroy the incoming missile is 
13%. An approximation of the natural single-pass efficiency of an FEL is only 
n=1/2N=2% so that a tapered undulator with high efficiency is needed. The tapered 
undulator described in the previous section would be more than adequate. This type of 


illustration is the essential driving factor in determining the FEL parameters for possible 


MW-class weapon systems. 
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B. LANL REGENERATIVE MW FEL AMPLIFIER 


The US Navy is investigating the possibility of using a MW-class Free Electron 
Laser for ship self-defense against anti-ship missiles (ASM). A design workshop resulted 
in two possible MW FELs, the oscillator and the regenerative amplifier. The MW 
oscillator design, although more compact overall, presents a challenge in recirculating an 
intense electron beam because of a phenomena known as coherent synchrotron radiation 
(CSR) feedback [11]. The regenerative amplifier design uses a longer undulator and 
relies on higher extraction efficiency to achieve high average power so that less current is 
needed from the accelerator and no recirculation is required. The regenerative amplifier 
design has been proposed by LANL and Boeing with the analysis being performed at 
NPS. This section explores the analysis of regenerative amplifier design and optimizes its 
efficiency. 


A simple schematic of the regenerative amplifier design 1s shown in Figure 3.1. 


____ Feedback Loop — 


Undulator 









Photo Injector 





Electron Dump 


Figure 3.1. Conceptual MW class Amplifier FEL. 


The amplifier design does not use energy recovery so there is no requirement to 
bend the high energy beam which can cause deterioration of the electron beam quality. 
However, without energy recovery, higher extraction efficiency is required. One- 
dimensional computer simulations are used to describe a single pass of the optical wave 


and to optimize the efficiency using a tapered undulator. The tapered undulator represents 
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an advantage when the electron phase acceleration 6 exceeds the deceleration that can be 
obtained without taper in strong optical fields. The single-pass FEL efficiency 7 of a 
tapered undulator is the fraction of the electron beam energy converted to laser light. The 
limit on 7 is determined by the maximum taper rate 6 that can maintain trapped 
electrons over the number of periods N in the undulator [6]. Gain G is the fractional 
change in power of the optical field in a single pass through the undulator. In weak fields 
without tapering, the FEL amplifier has high gain described by G(r) = (1/ gel 12)? V3¢] 

[12]. The gain is exponential in ¢ along the undulator with growth rate proportional to 


y°, where j=8N(ezKL) pF/y>mc’ is the dimensionless current density, 


p=3x10'1(A)/eczr; is the electron beam density, and F = "area of electron beam"/"area 
of light beam" 1s the filling factor. 


Table 3.1 shows the proposed parameters of the MW regenerative amplifier 


developed at LANL. 


Parameters Proposed Value 


Beam energy E 100 MeV 


Beam radius 7, 0.17 mm 


Energy spread Ay/y 0.02% 
Pulse duration 20 ps 
Pulse frequency 500 MHz 
Peak current / 400 A 
Average current / 0O.2A 
Undulator parameter K ey] 
Undulator length L 6m 
Number of undulator periods NV 300 
Undulator period 4, 2.0 cm 


Optical wavelength A 1.0 um 





Table 3.1. Parameters for the MW regenerative amplifier design developed at LANL 
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1. 1D Efficiency Optimization 


The proposed MW RAFEL uses a beam of 100 MeV electrons with a peak current 
of 400 A, yielding a peak electron beam power of 40 GW. The average current of the 
RAFEL is 0.2 A, so that the average electron beam power is 20 MW. The single-pass 
FEL efficiency 7 is the fraction of power extracted from electron beam in one pass 
through the undulator. The RAFEL will feedback on the order of 0.01% to 1% of the 
optical power and require an extraction efficiency of approximately 10% to 15% to 
provide 2 MW to 3 MW of optical power in the infrared (IR). 

In an attempt to achieve the desired efficiency, the initial undulator length is Z = 
4m corresponding to N = 200 periods. The dimensionless initial optical field a, and 
tapered undulator of strength 6 are varied to find the optimum efficiency. The maximum 
extraction efficiency found was 8% which is less than the nected efficiency. 

In an effort to improve the design efficiency, the length of the undulator is 
increased from 4 m to 6 m. The change in undulator length corresponds to dimensionless 


current density of j ~ 12,000 with N = 300 undulator periods. The values of a, and 6 


are varied to find the optimum efficiency using numerical simulations. The electron 


beam has a Gaussian spread o, =1.0, and the field strength is a, =150 with taper rate 
6 =2607 starting at time 7, =0.19 along the undulator. Figure 3.2 is the 3D plot of 


efficiency versus a, and 6 . The maximum efficiency 7,,,, 1s 13.3% with gain G = 556. 
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Efficiency Optimization 





Figure 3.2. 3D plot of Efficiency vs. Optical Field Amplitude 


and Tapered Field Strength. 
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The effect of changing a, away from its optimum value is smaller than changing 6 away 


from its optimum value as shown in Figure 3.3. 


Efficiency Optimization 


150} 
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Fig 3.3. Contour plot of Efficiency vs. Optical Field Amplitude 
and Tapered Field Strength. 
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For a better understanding of the physical mechanism underlying the operation of 
this MW RAFEL, examine Figure 3.3 closely in five regions. Region I is the "optimum 


region". As shown in Figure 3.4, the taper turns on at 7, = 0.19, and approximately half 
of electrons were trapped near resonance v, = 0, resulting in an acceptable efficiency of 
7 =13.3%. The final gain G,=556 corresponds to feedback of 
f =1MG, +1)=1.8x 10°. The large optical phase shift ¢ in the lower right-hand portion 


of Figure 3.4 represents good optical guiding which is an important part of a high gain, 
high current FEL. 





Figure 3.4. Phase space evolution for the MW RAFEL at optimum efficiency. 
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Region II is the "high tapering rate region". Increasing the tapering rate 6 beyond 
the optimum value, the efficiency drops to 8.3% due to fewer electrons being trapped near 
resonance v, =0 as shown in Figure 3.5. The taper turns on at 7, = 0.19 as before. The 
final gain G, =347 corresponds to feedback of f=1/(G, +1)=2.9x Loe g is 


relatively small which would cause less optical guiding. 





372 ) - 3n/2 0 c 1 


Figure 3.5. Phase space evolution for the MW RAFEL with high tapering rate. 
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Region III is the "low tapering rate region". For 6 below the optimum value, the 
efficiency drops to 6.5% due to insufficient taper rate as shown in Figure 3.6. While many 
electrons are trapped near resonance v,=0, the low taper rate does not extract 
significant energy from them. The taper again turns on at 7, =0.19. The final gain 


G, =294 corresponding to feedback of f=1/(G, +1)=3.4x 10°. 





Figure 3.6. Phase space evolution for the MW RAFEL with low tapering rate. 
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Region IV is the "high optical field region". Starting with high optical field 
strength a, , the separatrix does not have a chance to grow and trap electrons, resulting in 
a lower efficiency of 9.5% as shown in Figure 3.7. The taper turns on at 7, = 0.16. The 


final gain G, =144 corresponds to feedback of f=1G, +1)=6.9x10™. 





Figure 3.7. Phase space evolution for the MW RAFEL with high optical field strength. 
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Region V is the "low optical field region". Starting with low optical field strength 
a,, the initial separatrix 1s too small so that more electrons leak out from the separatrix, 
resulting in a lower efficiency of 10.9% as shown in Figure 3.8. The taper turns on at 
tT. =0.32. The final gain G,=25,700 corresponds to feedback of 
fAKGA AN) —3 9x0. 





Figure 3.8. Phase space evolution for the MW RAFEL with low optical field strength. 


The above results for efficiency should be considered an over-estimate since the effect of 


diffraction of the optical beam or longitudinal pulse effects are not included. 
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2. 3D Diffraction Optimization 


From Chapter II, the parabolic wave equation can be written as 


—] O a 
vie Z Jars, yzird=—il a) (3.1) 


Equation (3.1) governs the dynamics of the optical wave over many optical wavelengths in 
the longitudinal z-dimension, and the V*{ operator properly describes the diffraction of the 
optical wave in the transverse (x,y) directions. A coherent, freely-propagating optical 
wave of radius w, will spread due to natural diffraction as it traverses the undulator length 
L. The Rayleigh length, zw’/A, is the characteristic distance over which the initial 
optical mode area doubles in size. If the optical mode area becomes much larger than the 
electron beam, the optical field will decouple from the electron beam. Therefore, it is 
important that the Rayleigh length is comparable to the undulator length Z, and the optical 
mode waist w, is not too much larger than the electron beam radius r; [13]. 
Three-dimensional simulations are used to describe a single pass of the optical 
wave through the undulator. Initially, the optimal values from one-dimensional simulation 
were used as the inputs for a 3D simulation. The resulting 3D simulation is shown in 


Figure 3.9 with a, =150, 6 =260z, j =2.6x10", and the taper turning on at 7, = 0.2. 


The evolution of the optical mode, 





a(x, r) shown in the upper-left window of Figure 3.9, 
is dominated by diffraction after the taper is turned on at 7, =0.2. The intense electron 


beam provides guiding over the first one fifth of the undulator, but after tapering begins 
diffraction overcomes this focusing over the rest of the undulator. The top-center window 


shows a cross-sectional view of the optical mode, 





a(x, yl, at the end of the undulator, 


and this optical mode has a drop in optical power in the center. This drop is due to a loss 
of guiding of the optical mode [14]. The window in the top-right lists the dimensionless 


parameters for this simulation, o, and o, correspond to the radial size of the electron 


beam, while o, and o, account for beam quality. 
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Figure 3.9. 3D diffraction simulation. 


The graph at the right-middle shows the development of the optical phase along the 


undulator length. At 7, =0.2, the optical phase change along the undulator length 


becomes negative and drops rapidly which causes the optical wavefront to be excluded 
from the electron beam decreasing the interaction strength. The lower-right window is the 
plot of the dimensionless optical power growth along the undulator, and the plot of the 
natural logarithm of the single-pass gain along the undulator. At the middle-left shows the 
bunching current in the electron beam development along the undulator, with the end view 
in the center. The bottom-left plot is the electron beam phase velocity evolution along the 
undulator, and at the bottom center shows the electron phase space plot at the end of the 
undulator. Initially the optical guiding confines the optical mode until the taper turns on at 


t, =0.2, then diffraction effects begin to spread the beam. The optical phase begins to 


grow until the taper turns on, then the change in optical phase is negative as the guiding 
effect is lost. The bottom-center window shows that at the end of undulator, no electrons 


are trapped at resonance, and the tapering effect has failed. 
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In order to search for an optimum efficiency using numerical 3D simulations, the 
initial field a,, the taper rate 6, and the time to turn on the taper r, are varied. The 
optical field amplitude a, is varied from 50 to 275 in steps of 25. For each value of a,, the 
taper field rate 6 is varied from 30z to 3307 in steps of 257, and for each value of a, 


and 6, the turn-on time 7, is varied from 0.1 to 0.5 in steps of 0.02. First we explore the 


MW RAFEL efficiency by varying 6 and 1, at fixed values of ay. 


3D Efficiency Optimization (for a, = 50) 
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Figure 3.10. 3D and Contour plots of Efficiency vs. Tapered Field Strength 


and Taper turn-on time at a, = 50. 
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Figures 3.10 through 3.14 show the dependence of efficiency as 7, and 


6 are varied. These two parameters can have a significant impact on efficiency 
when moved away from their optimum values. The two graphs of each figure are 


plotted with the same data with the top one is in 3D and the bottom one is in 2D. 


3D Efficiency Optimization (for a, =100) 
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Figure 3.11. 3D and Contour plots of Efficiency vs. Tapered Field Strength 


and Taper turn-on time at a, = 100. 
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At the initial optical field amplitude of a, = SO, the efficiency peaks at 7. ~ 6.7% 
at values of 6 =2007 and +, =0.28. As a, increases, the efficiency peaks at slightly 
higher efficiency 77,,,. ¥ 7.0% at lower values of 6 =190z and +, =0.21 as shown in 


Figure 3.12. 


3D Efficiency Optimization (for a, =175) 
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Figure 3.12. 3D and Contour plots of Efficiency vs. Tapered Field Strength 


and Taper turn-on time at a, = 175. 
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The peak efficiencies 7,,,, continue to increase to a optimum efficiency 
Tmax ~ 1.2% Corresponding to optimum values a, = 225, 6 =180z, and taper turn-on 


time of 7, = 0.18 as shown in Figures 3.13. 


3D Efficiency Optimization (for a,=225) 
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Figure 3.13. 3D and Contour plots of Efficiency vs. Tapered Field Strength 


and Taper turn-on time at a, = 225. 
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Increasing optical field strength a, beyond its optimum value, the separatrix does 
not have a chance to grow and trap electrons, resulting in a slightly lower efficiency of 


7.1% as shown in Figure 3.14. 


3D Efficiency Optimization (for a,=250) 
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Figure 3.14. 3D and Contour plots of Efficiency vs. Tapered Field Strength 


and Taper turn-on time at a, = 250. 
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Next, explore the effect on efficiency by varying a, and r, at a fixed value of 6. 


As shown in Figure 3.15, starting at a low value of the taper field rate 5=55z, the 


maximum efficiency obtained was 3.2% with 7, ~0.2. This low efficiency is due to 


insufficient taper rate. 


3D Efficiency Optimization (for 6=557) 
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Figure 3.15. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 


and Taper turn-on time at 6 = 55z. 
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As 0 increases, the efficiency increases to 77,,,, ~ 5.3% at a slightly lower r, as 
shown in Figure 3.16. The peak efficiency increases to 7. ~ 7.2% corresponding to 


optimum values of 6 = 1807, a, =225, and an earlier taper turn-on time of 7, = 0.18. 


3D Efficiency Optimization (for 6=105r) 
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Figure 3.16. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 


and Taper turn-on time at 6 =105z. 


40 





As shown in Figure 3.17, when the optimum value of taper rate is used, 6 = 1807, 
varying the field a, has little effect on efficiency, whereas changing the taper turn-on time 


from its optimum value, r, ~ 18, results in a significant efficiency drop. 


3D Efficiency Optimization (for 6=1807) 
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Figure 3.17. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 


and Taper turn-on time at 0 =1807z. 
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Increasing 6 beyond its optimum value causes the efficiency to drop significantly 
due to fewer electron being trapped. As shown in Figure 3.18, the efficiency drops to 


only 1.9% at 7, =0.22 with 6 = 2807. 


3D Efficiency Optimization (for 6=2807) 
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Figure 3.18. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 


and Taper turn-on time at 6 = 2807. 
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Now consider the three-dimensional simulation of a single pass of the optical wave 
through the undulator at the optimum values of 6 =180z, 7, =0.18, and a, = 225. As 
shown in Figure 3.19, optical guiding occurs and stays with the intense electron beam 
through the undulator as seen in upper-left and middle-left windows. The top-center 
window shows a much tighter cross-sectional end view of the focused optical power than 
seen in Figure 3.9. The bottom-center window shows approximately 50% of the electrons 
are trapped, resulting in optimized gain and efficiency. The middle-right window shows 


the optical phase initially increasing until the taper turns on at 7, ~ 0.18, then decreasing 


slowly. The FEL performance in Figure 3.19 1s much improved over that of Figure 3.9. In 


steady-state operation, the final gain G,=616 corresponds to feedback of 


f =1G, +)) =16x10°. 
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Figure 3.19. 3D diffraction simulation at optimum values 
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IV. CONCLUSIONS 


The development of a speed-of-light hard-kill weapon system for military 
applications may represent a significant advancement in technology over present 
conventional kinetic weapon systems. The FEL has the potential for use as a weapon 
system against anti-ship and theatre ballistic missiles. The infinite magazine, rapid 
response, and wavelength tunability at high powers make it suitable and desirable for 
shipboard self-defense and for land-based theater ballistic missile defense. 

Three dimensional computer simulations indicate that the proposed LANL 


Regenerative MW FEL design does not yield the desired extracted efficiency of 7 ~15%. 
The optical field amplitude a,, taper rate 6, and taper turn-on time rt, were varied to 
search for the optimum efficiency, resulting in maximum efficiency of only 77, = 7.2%. 
Changing 6 or 7, away from their optimum values has a significant effects on the MW 


RAFEL design efficiency, but changes in a, do not have much of an effect. The efficiencies 
found here are still an over-estimate since they do not include longitudinal pulse effects. It 
remains to be determined whether a smaller efficiency 1s acceptable or if 15% efficiency is 
actually necessary in the FEL weapon design. If 15% efficiency is necessary, a longer 


undulator, perhaps 8 or 10 meters may be sufficient to reach that goal. 
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